Finite-size scaling exponents in the interacting boson model 
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We investigate the finite-size scaling exponents for the critical point at the shape phase transition 
from U(5) (spherical) to 0(6) (deformed 7- unstable) dynamical symmetries of the Interacting Boson 
Model, making use of the Holstein-Primakoff boson expansion and the continuous unitary transfor- 
mation technique. We compute exactly the leading order correction to the ground state energy, the 
gap, the expectation value of the d-boson number in the ground state and the E2 transition prob- 
ability from the ground state to the first excited state, and determine the corresponding finite-size 
scaling exponents. 

PACS numbers: 21.60.Fw,21.10.Re,75.40.Cx,73.43.Nq,05.10.Cc 



CO 



I 

o 
o 



(N 
> 
O 

(N 

m 
o 
in 
o 



-a 
c 

o 
o 



X 



The interest in the study of quantum phase transi- 
tions (QPT) has kept growing in the last years in dif- 
ferent branches of quantum many-body physics, rang- 
ing from macroscopic systems like quantum magnets, 
high-T^, superconductors Q] or dilute Bose and Fermi 
gases [2| to mesoscopic systems such as atomic nuclei or 
molecules Q. Although, strictly speaking QPT only oc- 
cur in macroscopic systems, there is a renewed interest in 
studying structural changes in finite-size systems where 
precursors of the transition are already observed £| . The 
understanding of the modifications on the characteristics 
of the QPT induced by finite-size effects is of crucial im- 
portance to extend the concept of phase transitions to 
finite systems. 

In the present study, we analyze these finite-size cor- 
rections in the Interacting Boson Model (IBM) of nuclei 
[l| , but the same technique can be applied to other boson 
systems, for instance to the molecular vibron model 
or to a multilevel boson model of Bose-Einstein conden- 
sates where similar QPT take place []]. 

The IBM is a two-level boson model that includes an 
angular momentum L = boson (scalar s-boson) and five 
angular momentum L = 2 bosons (quadrupole d-bosons) 
separated by an energy gap. The s and d-bosons repre- 
sent s and d wave idealized Cooper nucleon pairs. The 
algebraic structure of this model is governed by the U(6) 
group and the model has three dynamical symmetries in 
which the Hamiltonian, written in terms of the invari- 
ant (Casimir) operators of a nested chain of subgroups 
of U(6), is analytically solvable. The dynamical symme- 
tries are named by the first subgroup in the chain: U(5), 
SU(3) and 0(6). The classical or thermodynamic limit 
of the model was investigated by using an intrinsic state 
formalism that introduces the shape variables (3 and 7 
[lUGjl- Within this geometric picture the U(5), SU(3) 
and 0(6) dynamical symmetries correspond to spheri- 



cal, axially deformed and deformed 7-unstable shapes, 
respectively. Transition between two of these dynamical 
symmetry limits are described in terms of a Hamiltonian 
with a control parameter that mixes the Casimir oper- 
ators of the two dynamical symmetries. As a function 
of the control parameter, the system crosses smoothly 
a region of structural changes in the ground state wave 
function for finite number N of bosons. In the large N 
limit, the smooth crossover turns into a shar p O PT be- 
tween two well defined shape phases [9l llllll^ . ll3Lll^ . In 
particular, the transition from U(5) to 0(6) has been in- 
tensively studied in recent years because it has a unique 
second order QPT 0, Q, 0] associated with a triple 
point in the IBM parameter space [l4|. Furthermore, 
it was early recognized that the IBM Hamiltonian along 
this transition was fully integrable , and exactly solv- 
able [H El- 

Unfortunately it is difficult to use the exact solution 
to compute finite-size corrections analytically. Thus, we 
follow a different route which is based on the Continuous 
Unitary Transformations (CUTs) [H EH- Within 
this framework, we compute the first correction beyond 
the standard Random Phase Approximation (RPA) [2^ 
which already contains the key ingredients to analyze the 
critical point. As already observed in a similar context, 
[23I E3. l2fij . this 1/N expansion becomes, at this or- 
der, singular when approaching the critical region so that 
one gets nontrivial scaling exponents for the physical ob- 
servables (ground state energy, gap, occupation number, 
transition rates). In a second step, we take advantage of 
the exact solvability of the model to obtain numerical re- 
sults for large number of bosons which allows us to check 
our analytical predictions. 

Let us consider the U(5)-0(6) transitional Hamiltonian 
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Keeping terms of order (1/JV)° in the Hamiltonian ex- 
pressed in terms of the new 6's yields a quadratic Hamil- 
tonian which can be diagonalized via a Bogoliubov trans- 
formation. One then recovers RPA results [22I Ejf . 
At the next order (1/N) , the Hamiltonian is quartic 
and diagonalizing it clearly requires a more sophisticated 
method. To achieve this goal, we used the CUTs tech- 
nique 0, 120I l2l| . For an introduction to this method, 
we refer the reader to Refs. [23I I24I Ejf where CUTs were 
applied in a similar context. One introduces a running 
Hamiltonian 

H(l) = E (l) + A(l)n b + V(l):n 2 b :+W(l)P b i P b (3) 
+r(0 (p/ +P b )+ A(0 (P 6 + n 6 + n b P b ) , 

which is related to the initial Hamiltonian P(0) through 
a unitary transformation, namely H(l) — W (l)H(0)U(l). 
This transformation U is chosen such that H(oo) com- 
mutes with n b . In Eq. Q), : : denotes the normal 
ordered form of the operator O, and the notations for 
the 6's are the same as for the d's. The evolution of 
the running Hamiltonian is obtained from the flow equa- 
tion diH(l) = [r)(l),H(l)], where r?(/) = diW(l)U(l) is 
the anti-hermitian generator of the unitary transforma- 
tion. For the problem at hand, we consider the so-called 
quasi-particle conserving generator 29] , 



r,(l) = T(l) (Pi -P b )+ A(i) f P fe f n 6 - n b P ; 



(4) 



designed to ensure H(oo) commutes with n b , i. e. r(oo) 
A(oo) = 0. 



The flow equations can be solved exactly, order by or- 
der in 1/N, and the coefficients of the final Hamiltonian 
arc found to be 



where n d E M 4 d M ( with ^ = -S^ -1 ' ,!, 2 ) is the E (oc) 
d-boson number operator, Pj — s t 2 an d pt _ 

l^djjdl^ ; x is the control parameter that mixes 
the U(5) linear Casimir operator (x = 1) with the 0(6) 
quadratic Casimir operator (x = 0). The system un- 
dergoes a QPT at x c = 1/2, between a U(5) (spherical) 
phase for 1/2 < x < 1 and a 0(6) (deformed 7-unstablc) 
phase for < x < 1/2, when N -> 00 E3, Q II^I 

In the following, we shall restrict our analysis to the 
spherical (symmetric) phase which allows to investigate 
the critical point more simply than from the deformed 
phase. The Holstein-Primakoff boson expansion of one- 
body boson operators is especially well-suited to perform 
a 1/N expansion of the boson Hamiltonian Q. In the 
present case, it reads [26L |27| 
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where E(x) = x(2x — 1). One can then straightforwardly 
analyze the low-energy spectrum. 

The ground state of H(oo) is the state |0) with zero b- 
bosons, whose energy is Eq(oo). The first excited state is 
five-fold degenerate and corresponds to one quadrupole 
boson foL jO) , whose excitation energy is A(oo). These 
are the five components of the first 2 + excited state. 
For the two-boson states, things are a bit more compli- 
cated because of the W term which is not diagonal in 
( 2c ) the basis of states {6+ 6+ [0)} with n, v = -2, -1, 0, 1, 2. 



It is however easy to see that P b P b has a nontrivial ac- 
tion only in the subspace |&2^-2|0)' "Tf^o 1^) 1 ■ 
The corresponding 3x3 matrix has eigenvalues (twice) 
and 10. One thus finds that there are 14 degener- 
ate states with excitation energy 2[A(oo) + V(oo)], and 
one + state which is given by -^=Pj|0) with energy 
2[A(oo) + V(oo) + 5W(oo)]. Let us emphasize that the 
degeneracy is lifted at order (1/N) 1 , an effect missed at 
the RPA order. The 14 degenerate states are the 9 com- 
ponents of the first excited 4 + state and the 5 components 
of the second excited 2 + state. These 4+ and 2^ states 
are degenerate along the whole transition line because of 
the common 0(5) structure. Note also that, at fixed N, 
the state degenerates with the 4^~ and 2^~ in the U(5) 
limit. This low-energy spectrum is depicted in Fig. ^ for 
N = 40. The agreement between numerics and analytical 
results is pretty good and has been checked to improve 
when N gets bigger, as long as one is sufficiently far away 
from the critical point. Indeed, as can be seen in Q5IH[I . 
the 1/N order corrections diverge at x = 1/2. This sin- 
gular behavior already found in other models [23L I2T 
is a signature of the noninteger scaling exponents 
that we discuss below. 

The main strength of the CUTs is to allow the com- 
putation of expectation values of observables as well as 
transition amplitudes. Thus, one has to perform the 
unitary transformation of the observables in which one 
is interested. In the present case, all observables can 
be deduced from the knowledge of the flow of the op- 
erator &t (Z) = W(l)blU(l). For example, the average 
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FIG. 1: Comparison between analytical results (solid and dot- 
ted lines) and the numerical results (circle, triangle and square 
symbols) for the first excitation energies, with N = 40. 



number of d-bosons in the ground-state of the Hamil- 
tonian H is found as (rid) — (0| ^ M fejj, (c>o)6 At (oo) |0) . 
This quantity can also be computed using the Hellmann- 
Feynman theorem which yields (n c i) = -^[(1 + y)Eo] 
with y = x/(l — x). One then gets 



(nd) 



3x - 1 

2S(x)V2 

5x(l-x) 2 
16iV 



1 



(9) 



7x 



~(x) 2 E(xfl 2 



However, this theorem cannot be applied to compute 
nondiagonal matrix elements such as transition ampli- 
tudes. To illustrate the power of the CUTs for such a 
task, we focus on the B(E2) transition probability be- 
tween the ground state and the first excited state which 

in the standard 



is defined as B(E2) = Sl^OlQ^lO^)! 

J, M) basis, with = s^e? + dgS. The flow equations 
for b^l) can still be exactly integrated out order by order 
in 1/N and leads to: 
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The comparison between analytical and the numerical 
B(E2) transition probabilities for a system of N — 40 
bosons, is shown in Fig. [5] As for the excitation energies, 
there are divergences in the B(E2) values at the critical 
point, although they now appear even at the RPA order 
[28j. However, at finite N values no divergence should 
appear in the physical magnitudes or their derivatives 
with respect to the control parameter x, even at the crit- 
ical point. This obvious remark allows us to determine 
the nontrivial scal ing exponents. Such an analysis was 
proposed in Refs. [23J, |2J, |2j| , and we shall now briefly 
recall how it works. The 1/iV expansion of any physi- 
cal quantity $ has two contributions, the regular (reg) 



FIG. 2: Comparison between analytical results (solid dotted 
lines) and the numerical results (circles) for the B(E2) tran- 
sition probability, with TV = 40. 
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TABLE I: Scaling exponents for the ground state energy Eo, 
the gap A, the number of d-bosons in the ground- state (nd)GS 
and the B(E2) transition probability. 



and singular (sing) respectively, when x approaches the 
critical value x c — 1/2: 



(11) 



A close analysis of the singular part in the vicinity of the 
critical point x c shows that the singular part scales as 
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where T$> is a function depending on the scaling vari- 
able TVS (a;) 3 / 2 only. To compensate the singularity com- 
ing from E(x)^ <s> (or its derivative), one thus must have 
T^(x) ~ :r- 2 W3 so that $;»*(:,:„) ~ N -{n*+2&/3) _ In 

Table H] the computed scaling exponents for the low en- 
ergy physical quantities studied are summarized. 

In order to check these results, it is important to an- 
alyze the large N behavior of Therefore, we have 
numerically solved the problem by diagonalizing the bo- 
son Hamiltonian up to N = 1000. Details of this 
calculation will be given in a forthcoming publication 
[3l|. As can be seen in Fig. [31 an excellent agreement 
is found between the exponents predicted analytically 
and the numerical results. Let us underline that the 
scaling exponent for the ground state energy has been 
recently obtained by Rowe et al. [^3 by using the collec- 
tive model associated to the IBM Hamiltonian [22] . This 
mapping onto a quartic potential also explains why we 
found the same finite-size scaling exponent for the ground 
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FIG. 3: Plot of the singular parts of Eo, A, (rid) and 
B(E2)/(5N) at the critical point x c — 1/2, in a log 10 — log 10 
scale. 



state e nerg y and the gap (1/3) in other similar models 
[23L E3 . 123. However, such an approach does not allow 
to simply compute the finite N corrections and may not 
be suitable to obtain the behavior for observables such 
as B{E2). The CUTs method is thus, in this context, a 
very useful tool. 

In the present work, we have exactly computed finite- 
size corrections beyond the RPA in the symmetric phase 
of the IBM model. We have shown that the spectral 
properties at the critical point in the U(5)-0(6) tran- 
sition have well-defined asymptotic limits and we have 
calculated the TV-dependent scale factors. A natural ex- 
tension of this work would be to investigate the broken 
phase (x < 1/2) but the presence of Goldstone modes 
in the low-energy spectrum (at the RPA level) makes it 
more involved |3l| . 

The study of the scaling properties at the critical 
point of a QPT of a finite-TV particle model is of 
primer interest in several mesoscopic systems like 
nuclei, molecules, and other physical systems. The 
present results provide a tool to tackle such a study 
and to characterize the approach of the system to the 
critical regions as the number of particles goes to infinity. 
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